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The model of an approximate non-perturbative calculations in the SU(3) gauge theory is offered. 
This approach is based on the separation of initial degrees of freedom on ordered and disordered 
phases. The ordered phase is almost classical degrees of freedom, the disordered phase is completely 
quantum degrees of freedom. Using some approximations and simplifications for 2 and 4-points 
Green's functions an effective Lagrangian describing both phases from the SU(3) Lagrangian is 
obtained. The calculations show that ordered phase is squeezed by disordered phase into defects. 
These defects are: an infinite flux tube filled with longitudinal color electric and magnetic fields; a 
color electric hedgehog; a defect having either two color electric dipoles + two color magnetic dipoles 
or two color electric dipoles or two color magnetic dipoles. It assumed that the color defects are 
quantum excitations in a gauge condensate. The equations for the disordered phase are an analog 
of Ginzburg - Landau equation. 



I. INTRODUCTION 



In QCD exists an open problem of the description of non-perturbative degrees of freedom. In the most known 
kind this problem is connected to the problem of quark confinement. In this work we offer an approximate model of 
non-perturbative quantization. We separate SU(3) degrees of freedom on two parts: the first one are almost classical 
degrees of freedom (ordered phase), the second one are entirely quantum degrees of freedom (disordered phase or 
gauge condensate). The ordered phase is described using classical field equations for a small gauge subgroup. A gauge 
condensate is presented by scalar fields with corresponding scalar field equations. In order to obtain the equations 
describing disordered phase we use several assumptions and simplifications for 2 and 4-points Green's functions. 
Roughly speaking we suppose that 2 and 4-points Green's functions can be simplified by the following way: G2 ~ <t> 2 
and G4 as G\ + aiGi + ot2 (0(1,2 are constants). 

In the modern language we derive field equations for the dimensions two and four condensates. As a consequence 
we have the possibility to describe the gauge condensate on the microscopical level as a field distribution. Using this 
idea we can investigate color defects in the gauge condensate. In our investigation the color defects are some defects 
which destroy the global spatial homogeneity of the disordered phase. It can be color flux tubes, color electric dipoles, 
color magnetic currents and so on. 

The dimension two condensate \AV) has received a great deal of attention in the last few years, see for example 
Q-@ ■ The output of these investigations is that a non- vanishing condensate is favored as it lowers the vacuum energy. 
As a consequence of the existence of a non-vanishingcondensate (^4^), a dynamical mass parameter for the gluons 
can be generated in the gauge fixed Lagrangian, see Q , @ • 

An interesting idea is implemented in the so-called "cut-off" model G3 m which gluons having momenta 
smaller than a fixed value K are bound inside non-perturbative structures. This idea is similar to the one we will use 
in this work: probably these gluons are in disordered phase which we investigate here. 

In the other context the condensate notion arises in the dual superconductor model [l4| where the vacuum of a 
non-Abelian gauge model may be regarded as a medium of condensed Abelian monopoles. The monopole condensate 
squeezes the chromoelectric flux (coming from the quarks) into a flux tube due to the dual Meissner effect. In this 
work we show that the gauge condensate squeezes the color electric and magnetic fields into defects considered here. 

In Ref. [l5j it is shown that it is possible to arrive at an effective dual Abelian-Higgs model, the dual and relativistic 
version of Ginzburg-Landau model for superconductor, from SU(2) Yang-Mills theory based on the Faddeev-Niemi 
connection decomposition and the order-disorder assumptions for the gauge field. 
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II. HEISENBERG'S NON-PERTURBATIVE QUANTIZATION FOR QCD 

In this subsection we will apply a version of Heisenberg's non-perturbative quantization method 16] to QCD. The 
classical SU(3) Yang-Mills equations are 

q vT b^ _ Q ^ 

where T B V = d^A B - d v A B + gf BCD A%A B is the field strength; B,C,D = 1,...,8 are the SU(3) color indices; 
g is the coupling constant; f BCD are the structure constants for the SU(3) gauge group. In quantizing the system 
given in Eqs. Q - via Heisenberg's non-perturbative method one first replaces the classical fields by field operators 
A B — ■» A B . This yields the following differential equations for the operators 

d v T B ^ LV = 0. (2) 

These nonlinear equations for the field operators of the nonlinear quantum fields can be used to determine expectation 
values for the field operators A B , where (• • • ) = (Q\ ■ ■ ■ \Q) and \Q) is some quantum state. One can also use these 
equations to determine the expectation values of operators that are built up from the fundamental operators A B . For 

example, the "electric" field operator, £f = d A B — d z A B + gf BCD AQA B giving the expectation (£ B ). The simple 
gauge field expectation values, (A^(x)), are obtained by average Eq. (J2J over some quantum state \Q) 

Q 8 V T B ^ Q) = 0. (3) 



One problem in using these equations to obtain expectation values like (A B ), is that these equations involve not only 
powers or derivatives of (A 3 ) (i.e. terms like d a (A B ) or d a dp(A 3 )) but also contain terms like Q 3 ^ = (A B A^). 
Starting with Eq. (J3J) one can generate an operator differential equation for the product A B A^ thus allowing the 
determination of the Green's function Q B „ 

Q \A B (x)d y „T B ^(x)\ Q) = 0. (4) 

However this equation will in turn contain other, higher order Green's functions. Repeating these steps leads to an 
infinite set of equations connecting Green's functions of ever increasing order. This construction, leading to an infinite 
set of coupled, differential equations, does not have an exact, analytical solution and so must be handled using some 
approximation. 

Operators which are involved in equation are only well determined if there is a Hilbert space of quantum 
states. Thus we need to look into the question of the definition of the quantum states \Q) in the above construction. 
The resolution to this problem is as follows: There is an one-to-one correspondence between a given quantum state 
\Q) and the infinite set of quantum expectation values over any product of field operators, (x% , x% . . .) = 

(Q\A™(xi)A™(x2) ■ ■ ■ \Q). So if all the Green's functions - SS,"'." {%x , x% . . .) - are known then the quantum states, 

\Q) are known, i.e. the action of \Q) on any product of field operators AV? (xi)A™(xi) ... is known. The Green's 
functions are determined from the above, infinite set of equations (following Heisenberg's idea). 

Another problem associated with products of field operators like A™(x)Ay(x) which occur in Eq. © is that 
the two operators occur at the same point. For non-interacting field it is well known that such products have a 
singularity. In this paper we are considering interacting fields so it is not known if a singularity would arise for such 
products of operators evaluated at the same point. Heisenberg in his investigations of a quantized nonlinear spinor 
field repeatedly underscored that in a quantum field theory with strong interaction the singularities of propagators can 
be essentially smoothed out. Physically it is hypothesized that there are situations in interacting field theories where 
these singularities do not occur (e.g. for flux tubes in Abelian or non-Abelian theory quantities like the "electric" field 
inside the tube, (£%) < oo, and energy density e(x) = ((£f) 2 ) < oo are nonsingular) . Here we take as an assumption 
that such singularities do not occur. 



III. INFINITE COLOR FLUX TUBE DEFECT. 



In this section we will present a cylindrically symmetric color defect in a gauge condensate - infinite flux tube 
between dyon and antidyon. We reiterate that we have not found an exact analytical method for solving the full 
equations for all the Green's functions. Our basic approach in this case is to give some physically reasonable scheme 
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for cutting off the infinite set of equations for the Green's functions. In addition our approximate calculations will 
involve the decomposition of the initial gauge group into a smaller gauge group: SU(3) — ► SU(2) x £7(1). Physically 
we will find that this reduction splits the initial degrees of freedom (the SU (3) gauge potential components) into 
SU(2) x U(l) and coset components. The SU(2) x U(l) components will represent the ordered phase while the coset 
components, which will be represented via an effective scalar field, will represent the disordered phase. After these 
approximations one can perform analytical calculations which suggest the formation of a color electric flux tube. This 
is of interest since an important feature of the confinement phenomenon is the formation of electric flux tubes between 
the confined quarks. 

We suppose that quantum SU(3) gauge fields in quantum chromodynamics can be decomposed on a gauge conden- 
sate (ordered phase) and non-perturbative fluctuations (disordered phase). We suppose as well as that the flux tube 
is an axially symmetric defect in a gauge condensate. 



A. Reduction to a small subgroup 



First we define the reduction of SU(3) to SU(2) x U(l) following the conventions of Ref. [lOj. Starting with the 
SU(3) gauge group with generators T B , we define the SU{3) gauge fields, A^ = A B T B . SU{2) x U(l) is a subgroup 
of SU(3) and SU(3)/(SU(2) x U(l j) is a coset. Then the gauge field A^ can be decomposed as 



A 



A B T B _ A a T a 



A*T 8 + A™T m , 



A£ e 517(2), A* e 17(1) and A™ e SU(3)/ (SU(2) x U(l)) 



(5) 
(G) 



where the indices a,b,c... = 1,2,3 belongs to the subgroup SU(2) and m,n, ... = 4,5,6,7 to the coset 
SU(3)/ (SU(2) x U(l)); B are SU(3) indices. Based on this the field strength can be decomposed as 



*r-B rpB -r-a rpa 



-ITS rp» 



-77771 rpm 



where 



-77(1 




eSU(2), 




pa 
fits 


= d tl A a u -d v A;+ 9 r bc AlAl 


eSU(2), 




cf> a 


= gf amn A™A: 


eSU(2), 




J (JLV 


= h^+gf mn A™A: 


eU(i), 






= d^K - d v b^ 


eU(i), bp = A*, 






= F' m + G rn 


eSU(3)/ (SU(2) x 


u(i)) 


pm 


= <9,jA™ - d v A™ 


eSU(3)/(SU(2) x 


u(i)) 




= gf mnb (A^Al ~ A n v A\) + gf mnS (A^b u ~ A»6 M ) 


eSU(3)/{SU(2) x 


u(i)) 



where / are the structure constants of SU(2). The SU(3) Yang-Mills field equations can be decomposed as 



d v (F a ^ + $ a ^ l/ ) 



_ g famn A ra {pr^v + Qr^v^ 

-gf mnb [A™ (h b ^ + $ b ^) - a b v (F n ^ + G"^)] 
gf mnS \A n v W v - K {F n ^ + G"H] , 
d u (W +gf rnn A mf *A nu ) = -g Smn A™ (F n ^ + G"^) 



where d v [ 

j'rnnA J^n ^ 



• ] a = d v [- ■ ■ ] a + / abc a^[- • • ] c is the covariant derivative on the subgroup SU(2) and D v \ 
] A . The Heisenberg non-perturbative quantization procedure gives us 



D„ (f™^ + G miMV 



g jmnb 

gf mnS 



d u (h^ + 5/ 8mn A ro M m ') = -g 8 



'A 



A™ (h b ^ + 
Api^-bv (F nlMV 

~*~rn ^pnfj,v _|_ ^riftv 



(7) 

(8) 
(9) 
(10) 

(11) 
(12) 
(13) 
(14) 
(15) 

(16) 

(17) 
(18) 

(19) 

(20) 
(21) 
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B. Basic assumptions 

It is evident that an exact quantization is impossible for Eq's. H19|) - <|21|) . Thus we look for some simplification 
in order to obtain equations which can be analyzed. Our basic aim is to cut off the set of infinite equations us- 
ing some simplifying assumptions. For this purpose we will propose simple but physically reasonable ansatzen for 
the 2 and 4-points Green's functions - (A™{y)A^{x)), (^A a a (x)A b (y)A^(z)A^(u)S, (a% and 

(^A , ^(x)A , ^(y)AP l (z)Al(u)' K j - which are involved in averaged SU(3) Lagrangian. As was mentioned earlier we assume 
that there are two phases 

1. The gauge field components A a ^ 8 (A^ G SU(2),A* G U(l)) belonging to the small subgroup SU(2) x U(l) are 
in an ordered phase. Mathematically this means that 

<<£(*)) = K(x)) cl . (22) 

The subscript means that this is a classical field. Thus we are treating these components as effectively classical 
gauge fields in the first approximation. 

2. The gauge field components A™ (m = 4, 5,. ..,7 and A™ G SU(3)/ (SU(2) x U(l)j) belonging to the coset 
SU(3)/ (SU(2) x U(l)) are in a disordered phase (i.e. they form a gauge condensate), but have non-zero 
energy. In mathematical terms this means that 

(A™(x)) = 0. (23) 

According to lattice calculations the main contribution at calculations of any quantum quantities in QCD is 
brought with field configurations having Abelian monopoles. We use this fact to say that the contribution from 
the A™ component of gauge potential is negligible in comparing with the contribution from the A™ components 
M = 1,2,3 

(AZ(x)AZ(y)) « (A™(x)A?(y)) (24) 

where fi,i> = 1,2,3 are the spatial indices. After this remark we define an approximate expression for the 2-point 
Green's function as follows 

(A™{x)A n v {y)) = -V^f mpa f npb 4> a (x)4> b (y), (25) 
(A™{x)A n Q {y)) « (26) 

here </> a is a real SU(2) triplet scalar fields. Thus we have replaced the coset gauge fields in 2-points Green's 
function by effective scalar fields, which will be the scalar field in our effective SU(2) x fj(l)-scalar system. 

3. We suppose that the Green's functions with odd numbers of A™ are zero 

. odd number . 



^A™(x)...A:(y)J=Q (27) 

4. The correlations between the ordered (classical) and disordered (quantum) phases have the following forms 

(A™{x)A n v {x)b a {z)) = k 1 b a {z)(A^{x)A n v {x)), (28) 
(A^{x)AZ{y)b a {z)bp[u)) = klb a (z)b p (u)(A^(x)A:(x))+S mn M 1 .^b a (z)b (u), (29) 
(A^(x)A:(y)A a a (z)A b 3 (u)) = klA a a {z)A h p {u) (A™(x)A*(x)) + 6 mn M 2>flv A a a (z)Al(u), (30) 

where k\ t 2 are constants; M\^^ v are matrices which will destroy the gauge invariance of a final Lagrangian. 
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where \i,fj,\ are constants. 



5. The 4-points Green's functions can be decomposed by the following way. For A™ gauge potentials 
[A™{x)J%{y)Al{z)A*{u)) = X 1 (A™(x)A2(y)) ^(*)A«(u)) + , 

£ (5 mn v^ (A^(z)Al(u)) + 6™r, a0 {A™{x)A n v {y))) + 
^ (A:(y)Al(u)) + 8 n « Vvf} (A™(x)Al(y))) + 

£ (5 m ^ p (A n v (y)Al(u)) + 6 n P Vua (a™{x)A% (»))) +, 



t4 
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C. Derivation of an effective Lagrangian 



(31) 



For our quantization procedure we will take the expectation of the Lagrangian rather than for the equations of 
motions. Thus we will obtain an effective Lagrangian rather than approximate equations of motion. The Lagrangian 
we obtain from the original SU(3) pure gauge theory will turn out to be an effective SU(2) x U(l) Yang-Mills-Higgs 
system. The averaging SU(3) Lagrangian is 



where p,, v = 0, 1, 2, 3. According to item|2]of section ITlI Bl 



(32) 



(33) 



since = ~d^A + gf ABC A B A^ (we consider the stationary fields only), T A V = d^A v ~ d^A v + gf ABC A B A^ 
and correspondingly (T^ u T My ) has terms like (A B A m ) but (T A V T M ^) - (A B A B ^). ft means that the magnetic 
fields generating by the magnetic monopoles are more important by quantization in comparing with the electric fields. 
Therefore 



In Eq. (PH> J 70 "", and T rn ^ v are defined by equations 



(34) 



1. Calculation of (T^T^) 
We begin by calculating the first term on the r.h.s. of equation 1341) 



(35) 



Immediately we see that the first term on the r.h.s. of equation (|35|l is the SU(2) Lagrangian as we assume that A a 
and F® u arc effectively classical quantities and consequently 



/pa pafiu\ _ pa pupil* 



(36) 



The second term in equation (|35|l vanishes as F£ v is the antisymmetrical tensor but y&uv/ * s the symmetrical one 
(see Eq. (gSJ)). 

The last term which is quartic in the coset gauge fields can be calculated using (|31ll 



jamnppq ^ A n v A W A qV ) = -Ai ff 2 (0° ' (j) a - flf) 



Up to this point the SU(2) part of the Lagrangian is 



(37) 



(38) 
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2. Calculation of {T^ v T mtiv ) 

Next we work on the coset part of the Lagrangian 

{T™T m ^) = (F™F mfiU ) + 2gf mna <F™ (A au A n » - A a »A nv )) + 

2gf mnS (F™ (b v A n ^ - b»A nu )) + g 2 f mna f mpb ((A a „A n ^ - A^A") (A bp A p " - A bfl A pu )) + 
g 2 f mnS f mpS ((bvA^ - b^AZ) (b u A p " - b p A p,J )) + 

g 2j m najmpS - A«A™) (b U A P » - b^ A pU )) . 



(T™T m ^) = - 4 (d^ a ) W) ~ 2gk 2 f abc (d^ a ) A b ^ c - g 2 k 2 - 3g 2 k 2 {<j> a <j> a ) 



(39) 



After the calculations we have 

(F™F m ^) = -4(d^ a )(d p r), (40) 

g 2 f mna f mpb ({A a v Al - A*A$) {A bv A w - A bp A pv )) = -g 2 k\ A^A a " + 2 (M 2 )^ A a ^A% (41) 

g 2 fm n 8fmp8 ^ bvA n „ 6 ^«) ^ ^ _ ^ = -3 ff 2 fc 2 (</>> Q ) b^ + 2 (m 2 )^ (42) 

2gf mna (F™ (A av A n »- A afi A nu )) = -2gk 2 f abc {d^) A bp <j) c , (43) 

g 2 f mna f mpS ((A^- A^Afj (b"A Pfi -b»A pi ')) = 0, (44) 

2gf rnn8 (F™ (b u A np -b p A nu )) = (45) 

where 

(M 2 )^ = 9 2 W-Mf), (46) 

(m 2 )^ = 3 5 2 {M^ a rf v - Mf) . (47) 

In what follows we assign these matrices as diagonal 

(M 2 Y V = M^ v , (48) 

( m y = m^f". (49) 

on fj, there is no summing, M 2 and m 2 are constants. Finally 



2(m 2 r^ + 2(Af 2 r^. 



(50) 



3. Calculation of (T^T 8 "") 

We end by calculating the second term on the r.h.s. of equation l|34|) 

(Fl v F 8pv ) = (hpuM*") + 2g {h^A m ^A nv ) + g 2 f™fp* (A™ A n v A pp A qv ) . (51) 

The first term on the r.h.s. of l|51|l has the classical field consequently 

(V^H = K u h^. (52) 

The second term in equation (|51|) vanishes as is the antisymmetrical tensor but (^A m ^A nu ) is the symmetrical one 
(see Eq. (|2SJ)). The last term which is quartic in the coset gauge fields can be calculated using lj3T)l 

9 2 / 8 "/ 8M (A^A^A^A^) = ®-\ l9 2 - • (53) 

Thus 

(FlF 8 ^) = + h l9 2 - fif ■ (54) 
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D. An effective Lagrangian 

Thus after the quantization A™ degrees of freedom we have 

{F£ V T A * U ) =F^ v F a ^ v + h^h^ - 4 (d M 4f l ff i 4> a + ^ge^d^A^cf) + i^Jfe|A* A 6 *V c c ) + 

V 2 4 J (55) 

\\ l9 2 -AY- 3kW (O a ) b^V + 2 (Af 2 r + 2 (m 2 )^ b,K 

After the redefinitions </> a — > ^= , /ii — > ^= , 9g 4 Al — > Ai we have the full averaged Lagrangian 

£ eff = -\ (T%T A » V ) = - \f%F^ - \h^ v + \ (D,r) (><r) ~ l -g 2 k 2 2 AlA^<t> b ~ 

£ {rr — Mi) + 6 -kW b^v - \ (M 2 r a*a* u - \ ( m 2 r 

where D^<f> a = d^cf) a + ^k2e abc A bti <fi c . The original pure SU(3) gauge theory has been transformed into an SU(2) x U(l) 
gauge theory with broken gauge symmetry and coupled to an effective triplet scalar field. It is necessary to note that 
the coupling constant for an effective triplet scalar field is ffo not g. This Lagrangian we will investigate numerically 
on the next sections. 

Our final result given in l|56[) depends on several crucial assumptions. In the next section we will investigate a 
non-topological flux tube solution to this effective Lagrangian. 



E. Initial equations 

In this section we will investigate a non-topological flux tube solution to the effective Lagrangian (|56|l . i.e. a 
cylindrically symmetric defect in a gauge condensate. We will see that the corresponding solution is an infinite flux 
tube filled with longitudinal and transversal electric and magnetic fields. It allows us to say that the obtained solution 
is the flux tube stretched between two color charge with chromoelectric and chromomagnetic fields. The separation 
between charge is infinite one. 

We will start from the SU(2) x U(l) Yang - Mills - Higgs field equations of the Lagrangian i|56[l with broken gauge 
symmetry 

D v F allv = ge abc (l) b D^(f, c - (M 2 )^ A a v - ^k 2 g 2 A b ^(f, b (j) a , (57) 

d v h» v = |*?</ a (m 2 Y v b Ul (58) 

= -\0 a - vl) + Z -g 2 k\b^^ - X -g 2 k\A\A a ^ b (59) 

here F afiy = d^A® - d v A a ^ + g] abc A b ^A c v is the field tensor for the SU(2) gauge potential A£; a,b,c= 1, 2, 3 are the 
color indices; D v [- ■ ■ ] a = d v [- ■ ■ ] a + gf abc A b ^[- • • ] c is the gauge derivative; <j) a is the Higgs field; Ai, g and Pl are some 
constants; (M 2 ) M!y and (ro 2 )^ are masses matrices which destroys the gauge invariance of the Yang - Mills - Higgs 
theory. In the present case we are postulating that the masses of the two gauge bosons are different. This was done 
in order to find electric flux tube solutions for the system l|57|) - <|59|) . For equal masses electric flux tubes did not exist 
for the ansatz used. In numerically solving Ij57|l - (|59|) the two mass come out close to one another, but are not equal. 
We will use the following ansatz 



9 9 9 



&„ = ^{h(p),0,0,w(p)}; (M 2 )^ = diag{Mo,Mi,0,0}; (m 2 )^ = diag |m , 0, 0, ^ 
here z,p,ip are cylindrical coordinates. The color electric and magnetic fields are 



(60) 



El = -fv, El = -If, H% = -v', (61) 
9 9 9 

El - hf, HI - 1* (62, 
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Substituting this into the Yang - Mills - Higgs equations (|57[1 - <|59|) gives us 



here we redefined <p/(po 
WO,3/0O -* ™o,3, P<Po -> x; <pQ = </>(0) 



// 

V 


V 
X 






f" 


£1 

+ L = 

X 


/( 


4> 2 + v 2 




<t>' 

+ — = 

X 


<p 


-f + ' 


h" 


ti 

+ — = 

X 




[hu 2 


w" 


w' 

X 


w 




* 0, 


f/<h - 


■+ /, v/(f> 



Mi) 



kf -h z 



/, v/(f> -> u, /i/(/> -> /i, tu/^o -> Mi/0o -> Mi: M 0il /(p 



(63) 
(64) 
(65) 
(66) 
(67) 
M ,i, 



F. Numerical investigation 



For the numerical calculations we choose the following parameters values 



0(0) = 1, Ai = 0.1, u(0) = 0.5, 



(68) 



We apply the methods of step by step approximation for finding of numerical solutions (the details of similar calcu- 
lations can be found in Ref. [ll|). 

Step 1. On the first step we solve Eq. ijfiSl) (having zero approximations fo(x) = — ° 2 2 / ^ ) , 4>o{x) — 1.3 ° 2 3 / - r . 

cosn ^ J cosli y — J 

The regular solution exists for a special value Mj" i only. For Mi < Ml . the function Vi (x) — > +oo and for Mi > M* 4 
the function i>j(x) — > — oo (here the index i is the approximation number). One can say that in this case we solve a 
non-linear eigenvalue problem: v*(x) is the eigenstate and M*^ is the eigenvalue on this Step. 

Step 2. On the second step we solve Eq. Qoty using zero approximation (f>o{x) for the function f(x) and the 
first approximation v1(x) for the function v(x) from the Step 1. For Mq < Mq 1 the function fi(x) — * +oo and for 



M > M * j the function fi (x) 



Again we have a non-linear eigenvalue problem for the function fx(x) and 



o.i- 



Step 3. On the third step we solve Eq. I|65|l using the first approximations f*(x) and v\(x) from the Steps 1,2. 
For fi\ < fil 1 the function cj>i(x) — > +oo and for [i\ > 1 the function 4>i(x) oscillates and tends to zero. Again we 
have a non-linear eigenvalue problem for the function <pi(x) and (Iq 1 . 

Step 4. On the fourth step we repeat the first three steps that to have the good convergent sequence 
Vi(x), fi(x), 4>*(x). Practically we have made three approximations. 

Step 5. On the fifth step we solve Eq. I|66l) using the first approximations 4>*(x) from the Step 3. This equation 
is exactly the Schrodinger equation for the function h(x) with the potential 4> 2 (x) and eigenvalue to 2 ,. 

Step 6. On the sixth step we solve Eq. I|67|l using the first approximations <fi*(x) from the Step 3. For 7773 < nij 1 
the function wi(x) — > +00 and for 7773 > 777,3 1 the function wx(x) — > —00. Again we have a non-linear eigenvalue 
problem for the function Wi(x) and 777.3 i- 

Step 7. On this step we repeat Steps 1-6 necessary number of times that to have the necessary accuracy of definition 
of the functions v*(x), f*(x),(j>*(x),h*(x), w*(x). 

After Step 7 we have the solution presented in Fig's. H - EU These numerical calculations give us the eigen- 



values M ( 

«*(*), ro 



re 1.2325025678, Ml 

(f>*(x),h*(x),w*{x). 



1.18060857823, 



1.18356622792, 



1.29306895 and eigenstates 
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FIG. 1: The functions /* (a;), v*(x), h*(x), w* (x), (j)*( x ) 




FIG. 2: The longitudinal chromoelectric E^(x) and chromomagnetic H^(x) fields. 



It is easy to see that the asymptotical behavior of the regular solution of equations (|63|) - (|67|1 is 

exp j-xy / 2A^uf} 



<P{x) ~ Mi 



cx P {-xVmi-^o} 
j ■ 

expj-xVMi - M i] 



exp 



{-a^/i? - rug} 



exp |— x\/ pf — 



where 0oo, Wqo, /qq, /loo, Wo, are some constants. 

Let us note the following specialities of the solution: 

1. The flux of the chromoelectric field is nonzero 



$ £ = 2?r / E%pdp = 2tt / fvpdp ± 0. 
o o 



(69) 
(70) 

(71) 
(72) 
(73) 



(74) 
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-0,14 J 



FIG. 3: The transversal chromoelectric Ep :8 (x) and chromomagnetic H^{x) fields. 



2. The flux of the chromomagnetic field is zero 

oo oo 

<&H = 2tt J H*pdp = 2tt J ^-dp = 2tt [to(oo) - w(0)] = 0. (75) 



3. The longitudinal electric field 

is essentially non-Abelian field. 

4. The longitudinal magnetic field 

is essentially Abelian field. 



El = gA}Al (76) 



p dp 



5. Comparing the obtained flux tube with the Nielsen - Olesen flux tube we see that the scalar field <fi a is an analog 
of the Ginzburg - Landau wave function for Cooper pairs and Eq. (|59|l in QCD is the analog of the Ginzburg - 
Landau equation. 

6. According to the previous remark the scalar field <j) a probably describes a monopole condensate. 



IV. A DEFECT WITH TWO INFINITESIMALLY CLOSED ELECTRIC AND MAGNETIC DIPOLES. 

In the previous section we have obtained the color defect (infinite flux tube with the longitudinal electric field) 
where the color charges are on infinite distance from each other. In this section we want to investigate the case with 
infinitesimally closed color charges. 



A. SU(3) — > U(l) + coset decomposition 

Now we will consider the case when A°' m (a = 1, 2, 3, m = 4, 5, 6, 7) degrees of freedom are quantized but remains 
in a classical phase. Again we start from the initial SU(3) Lagrangian £ su(3) where we have to average over A®' m . 
In order to do this we have to add assumptions in addition to the section UllBI 

1. On this step the gauge field components A® (A^ S SU(2)) are in a disordered phase (i.e. they form the gauge 
condensate), but have non-zero energy. In mathematical terms this means that 



(a;(x)) = o. 



(78) 
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The reasonings similar to those which were resulted regarding section UlI Bl result us in expression 

(a^aKv)) = ~i hv r m f pn r{xW(y) (79) 

here (f> m ,m — 4,5,6,7 are real scalar fields. Thus we have replaced the SU(2) gauge fields by effective scalar 
fields, which will be the scalar fields in our effective E/(l)+(two scalar field) system. 

2. The 4-points Green's function for A a gauge potential can be decomposed by the following way 



(A^x)Al(y)A c a (z)A$(u)) 



\ 2 y{A1{x)A b M) (A c a (z)Af,(u)) + , 

£ (*»V (Al{z)A$(u)) + S^naf, (A«(x)Al(y))) + 
^ (6 ac Vfia {A b v {y)A${u)) + S bd Vv0 (A;(x)A c a (y))) + 
^ (6 a %p (A b u (y)Al(u)) + 6 bc Vl/ct {Al{x)A${y))) +, 
y| {6 ab V»J cd V af 3 + S ac V m S bd V.p + 5 ad r!^5 bc Vva ) 



(80) 



where A2,/ia are constants. 

3. The correlations between the ordered (A & = b^) and disordered (A®' m ) phases have the following forms 

(A^ m (x)A b / l (x)b a (z)) = k 1 b a (z)(A^ m (x)A b ^(x)), (81) 
{A^ m {x)AY{y)b a (z)b {u)) = klb a (z)b p (u) (A;'"\x)A b j n (x)) + S^M^b^b^ul (82) 

where k% are constant and M\ iliv are matrices. 

4. We suppose that the Green's functions with odd numbers of At are zero 

odd number 



^(x)...A b u (y)J=0 (83) 

5. The correlations between A't and A™ quantum degrees of freedom are 

{A™{x)A n v (y)A a a {z)A b (u)) = kl{A™{x)A n v {x)) {A%{z)A h p {u)) . (84) 
Thus we have to calculate the following effective Lagrangian 

C eff « -\ (T* V T A » V ) = ~\ ((^ a n + {mT m ^) + (^ 8 ^» (85) 
taking into account that A c ^ ra are in disordered phase but A^ remains in ordered phase. 



B. Calculation of {T* v T atlv ) 



The calculations are similar with the section UlI G II We begin by calculating the first term on the r.h.s. of equation 

(T^T a ^) = {{d^A a u - d v Alf) + 2gf abc ({d^Al - d u A«) A b »A c ») + g ^f^ f ade (^^^ + 
2gf amn (F a A m ^A nv ) + g 2 f amn f a Pi (^A m A n A vtl A qv ) 

The first term is 

{{d^Al - d„Alf) = -Id^d^iT (87) 
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The second term ((d^AI — d„A°i) A hil A cv ^ on r.h.s. of equation (|35|l vanishes as it has odd number of A" The term 

(F^A™ 11 A nu ) is vanishes as F£ v is the antisymmetrical tensor but (A mf2 A nl/ ) is the symmetrical one (see Eq. (JZHl)- 
The two terms which are quartic in the coset gauge fields can be calculated using lj3*T)) and (|SU|) 

g 2 f a bcf ade ^b^^ = ^ g 2 _ fif ^ ( g g) 

famnjapq ^ A n v A w A qv ) = ^X ig 2 ((j) a (j) a - /I 2 ) 2 . (89) 

Up to this point the first term {T^T 11 ^) of the Lagrangian i|85|) is 

(j%s" tv ) = -\d^ m d^ m + h l9 2 (rr - M?) 2 + \w (<r<r - m 2 .) 2 . (90) 



C. Calculation of {J^JT 



The next term is 

(T^T™^) =({dpA? - d u A™) (d^A mv - d v A m ")) + 2gf mna ((d^A™ - d v A™) (A aL, A nf * - A a "A nv )) 



2gf mn8 (OVl™ - d v A™) {b v A n » - b"A nu )) 



g 2 fm na fm pb ^ ( A a A n _ j^aj^ (jbv ^ ) + ( 91 ) 

g 2 f mn8 f mp8 {{b v Al - b^Al) {b u A p » - WA pv )) + 
5 2 / mna / mp8 {{K A l ~ ip VAW - b fJ "A pu )) 

After the calculations we have 



((d ll A^-d v A^)(& i A mv -ff , A m i i )) = -A(d^ a ){d^ a ), (92) 

gljmnajmpb ((^°^» _ A^A™) (A bv A p " - A bfl A pu ) ) = ^g 2 k 2 2 ((f) a <j) a ) (</>"> m ) , (93) 

g 2 f mnS f mpS _ b ^ ^ ^ _ = ^2^ (0>Q) ^ + 2 ^ (g4) 

2gf mna ((d^A™ -d u A™) (A av A n » -A a »A nu )) = 0, (95) 

g 2jmnajmpS ((A%A% - A^A") {b" A Pt * - b^ A P ")) = 0, (96) 

2gf mnS ((d^A™ - d v A™) (b v A n "-b fl A nv )) = 0. (97) 



Finally 

(Jftj™"') = - 4 (d^ a ) W) + \g 2 k 2 2 (rr) {<f> m 4> m ) - 3g 2 k 2 (<f> a r) + 2 (m 2 )^ b,K (98) 



D. Calculation of (•7 r *„.7 r8 '" / ) 

We end by calculating the third term on the r.h.s. of equation l|85|l 

{T^T 8 ^) = ( + 2 3 {h„ v A m "A nv ) + g ^f mn fpi (A™ A n v A w A qv ) . (99) 
The first term on the r.h.s. of H51|l has the classical field consequently 

(V/O - K^ V - ( 10 °) 

The second term in equation (j51(l vanishes as h^u is the antisymmetrical tensor but ^A m ^A nu ) is the symmetrical one 
(see Eq. I|25|l). The last term which is quartic in the coset gauge fields can be calculated using j5J 

g 2j8ranjS Pq (A™A™A P »A qV ) = \\ X g 2 (cj) a (j) a - ^ . (101) 



Thus 



(J*,?*"") = h^hf" + ^A l5 2 - Mi)' • (102) 
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E. An effective Lagrangian and field equations 



Thus after the quantization A®'" 1 degrees of freedom we have 
-4 (C eff ) = Vft"" " \ (d^ m ) (d^ m ) ~ 4 (d^ a ) (W) + — Mi) 2 + J5 2 A 2 (<T<T - Ms) 2 + 

9 - 9 2 ki (rr) {<r<T) - 3.9 2 fc 2 ^) b ^ + 2 (m 2 )^ 



(103) 



After the redefinition 6 m — » -^6 m , d) a 

v3 



^f,M2 -> ^jM2,Mi ^Mi,4g A 2 -> A 2 , igff A x -> A x we have 

1 
2 



^1 
4 



2 _ A 
"4 



■ 4> ~M 2 ) 



The field equations for the U(l)+(two scalar) system are 



8»8 U 



•Ah 



2 „2 



+ A 2 (<T<F - /i 2 ) 



d v h^ = ^ 2 fc 2 (^ a )6"-(m 2 )^6, 



(104) 



(105) 
(106) 
(107) 



F. Numerical solution 



We will use the following ansatz 

4 




{f(r,9),0,0,v(r,9)}; <j> a 



diag <^ mo, 0,0, 



1 2 

9 3h 



<f>(r,6); <f 



1 1 



9 k 2 V3 



x(r,ey, 



r 2 sin 2 1 



here r, 9, ip are spherical coordinates. The corresponding color (but Abelian) electric and magnetic fields are 

= if^i H s = 1 dh * 
8 r d6' 9 r sin 9 dr ' 

db t „ 8 1 db d 



E 



8 



dr 



'■ sin 6* 89 



Substituting (|108[) into the U(l) - Higgs equations p05|) - (|107|l gives us 



1 d ( 2 90 



x- 



2 <9x \ dx 



1 9 

sin 

z 2 sm0 9(9l "89 



X 2 + Ai 



1 9 f sin ^ 
x 2 dx V dx ) ' x 2 sin 9 89 \ 89 



1 8 { 2 8 X 

x 



(0 2 " Mi) - (/ 2 
X [0 2 + A 2 (x 2 ~ Ml)] , 



x 2 sin 2 i 



i a / 2 a/\ | i 8 ( sin0 9f . = , 

x 2 <9x \ 9x / x 2 sin 9 89 \ 89 



8 2 v sin 9 8 ( 1 8v\ 
'dx 2 + ~89 \^n~989J 



"/V 


)% 2 - ? n 2 








2 




1 </> 2 -™ 2 



(108) 



(109) 
(110) 

(111) 
(112) 

(113) 
(114) 



here we redefined 0/0 o -> 0, x/0o -> X, //0o -> /, «/0o -> u, h^Mi -> Mi, ^%M2 -> M2, ^tAi -> Ai, 
^|p-A 2 — > A 2 , r0o — > </>o = 0(0)- The equations (|lll|l - (|114|l are partial differential equations and it is not quite 
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clearly how they can be solved numerically to find eigenvalues (11.2 and 777,0,3. I n order to avoid this problem we 
assume that the effect from the U(l) Abelian field is small and then 



f(x,6) = f(x) (cos0-^cos 3 
v(x,9) = v(x) sin 2 9 cos 9. 
In addition we have to average Eq. over the angle 9 using the following average 



- [ f 2 (x, 9) sin 9 d9 = ^-f(x) 7 

7T J bill 

4 v 2 (x) 





7T 

1 f v 2 (x 

7T 



■sin9 d9 = 



(115) 
(116) 
(117) 

(118) 
(119) 



157T X 2 

Finally we have the following ordinary differential equations which can be numerically solved as an eigenvalue problem 



X 

2 , 

-x 

x 
12 



X 2 + Ax U 2 - M ?) 



( 



f 2 -—- 

\^637T 157T X 2 



X W + A2 (x 2 - l4)] , 
2 



/" + -/' - sf = f 



-) -m 



— 777 



For this choice we should have 

Ee(x,9)\ e=0 „ = 0, H e (x,( 

E r (x,6)\ x=0 = 0, H r (x,( 

The series expansions near x — for the functions are 



= 0. 



10 = 0, 7T 









x 2 


<P(x) 


= 00 + 


03 


2 








x 2 
2 


X{x) 


= Xo + 


X3 


m 








= /3-jr 
6 


+ 




v(x) 


X 3 






= v 3 - 


+ 





(120) 
(121) 
(122) 

(123) 



(124) 
(125) 



(126) 
(127) 
(128) 
(129) 



provide the constraints (| 124(1 (|125|l . We will search for a regular solution with the following boundary conditions: 



0(0) = 1, 0(00) = m, 

x(o) = xo, x(°°) = 0, 

/(0) = /(oo) = 0, 

v(0) = v{oo) = 0. 



(130) 
(131) 
(132) 
(133) 



The numerical calculations are carried out as well as in the section UlI Fl with rj- = 6. The result is presented in Fig 
Hand /ji = 1.617525 . . . ,/x 2 = 1.4930404 . . . , m 3 = 3.46539926281310275 . . .. 
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FIG. 4: The functions f*(x),v*(x), X*0)> <j>*(x) 



It is easy to see that the asymptotical behavior of the regular solution of equations (|120fl - l|123|) ^ s 

exp ^-x y / 2X 1 ^ 



(x) PS /Xi 



exp {-xa//^ - A 2 ^| 



exp < —a; 



(fe) Mr 



u(x) w Woo exp ^ -x^/ ( ) /xf-mj 



where 0oo , Xoo , /oo , ^oo are some constants. 

In Fig. (JSJl and JBJ the distribution of color electric and magnetic fields are presented. 




(134) 
(135) 

(136) 
(137) 




FIG. 5: The profile of the Abelian electric field E s 



FIG. 6: The profile of the Abelian magnetic field H s 



V. COLOR ELECTRIC HEDGEHOG 



In the previous section the defect in the gauge condensate is considered which represents two infinitesimally closed 
color electric dipoles and two color magnetic currents. In this section we want to investigate an isolated electric 



16 



hedgehog in the gauge condensate. It is not pure electric charge as at the origin the radial electric field is zero. 
The origin of the electric hedgehog is the nonlinear interaction between ordered and disordered phases (nonlinear 
interaction between gauge potential components in the initial SU(3) Lagrangian). 

The initial equations describing such defect are Eq's (|lllfl - l|114[l with v — and we consider spherically symmetric 
case / = f(r), 4> = 4>{r), \ = x(r) 



x 

x" + -x' 

x 

/"+-/ 

X 



[X 2 + Ai (0 2 - M?) - f 
X [4? + A 2 ( X 2 - i4)] , 



(138) 
(139) 

(140) 



The numerical calculations are the same as in the previous sections and the profiles of gauge components and 
radial electric field are presented in Fig's [7| |H1 In this case /ii = 1.6098972 ... ,M2 = 1.4854892 too = 
2.9355398046322565 ... 
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/^T 6 


3 10 
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E r (x) 




-0,10- 
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FIG. 7: The functions (p*(x), X*{x), f*(x) 



FIG. 8: The radial electric field E r 



VI. TWO COLOR ELECTRIC DIPOLES 

In this section we want to investigate Eq's (|lll|) - i|114fl with v = 0, i. e. without magnetic field 

Mi) - f] , 



1 


d 




X 2 


dx 


1 






X 2 


— 

dx 


1 


d 


\ dx 


X 2 


dx 





1 


d 


X 2 


sin t 


)d6 




1 


d 


X 2 


sin £ 


we 




1 


d 


X 2 


sin£ 


)d9 





m 




x 2 + Ai (4? 


sin ( 


do ) 










<^ 2 + A 2 (x 2 


sin( 


de ) 


= x 




m 


= f 




sin( 


de J 


f(t)> 



f(x,e)=f{x) cos 0-- cos 3 



and 



Again we assume that the effect from the U(l) Abelian field is small and then 

4>{x, 9) « 4>{x), x{%, 8) ~ x( x )- 



(141) 
(142) 

(143) 

(144) 
(145) 
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Thus we will numerically investigate the following equations set 

2 



6" + -</>' 

x 



X +Ai 



-X' = X[0 2 + A 2 (x 2 -/^)], 



12 



/" +~f'--of = f 



(146) 
(147) 
(148) 



The result is presented in Fig. Hand m = 1.616718 . . . , M2 = 1.4926212 . . . , m = 3.670495292 . . . The field distribution 
of Abelian electric field E s in this case is the same as in Fig. [5] and similar to the field distribution between two 
electric dipoles lying on the line. 




FIG. 9: The functions <j>*{x), x* (x), f* (x) for two color electric dipoles. 



VII. TWO COLOR MAGNETIC DIPOLES 



In this section we want to investigate Eq's I|lll[l - qil4jl with / = 



1 9 f 2 d£\ , 1 d 
x 2 dx \ dx ) x 2 sin 9 d8 

1 d ( 2 d x \ 1 d 
x 



sm( 



89 
,dx 



x 2 + M 



Mi) 



x 2 sin 2 t 



x 2 dx \ dx J x 2 sin 9 89 \ 89 
8 2 v s'm9 d ( 1 dv 



sin^ = XW + M{X 2 -1%)], 



dx 2 



89 \ sin 9 86 



3fei 
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and 



v(x, 9) = v{x) sin 2 9 cos 9. 



(149) 
(150) 

(151) 
(152) 



Doing the same assumptions about smallness of the U(l) Abelian field we have again the approximations l|145|) 
and the following equations set 



4 v 2 

157T X 2 



X" + -X' = X[0 2 +A 2 (x 2 -^)], 
x 



-V = V 



J.2 2 

<P - m 3 



(153) 
(154) 
(155) 
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The result is presented in Fig. [TOJ and m = 1.6179587 ...,fi 2 = 1.4931567 . . . , mi = 3.465550906326 . . .. The field 
distribution of Abelian magnetic field H 8 in this case is the same as in Fig. EJand similar to the field distribution 
between two magnetic dipoles. 




FIG. 10: The functions cj>*(x),x*(x),v*(x) for two color magnetic moments. 



VIII. PARTON-LIKE DEFECTS IN A GAUGE CONDENSATE 



In this paper we would like to include quark field. The initial effective Lagrangian is 



2\ 2 



\g 2 k 2 (rr) (<r <n + \g 2 k\ (rr) ~ x - (m 2 y v + $ [r + <?a« a 8 ) - m j v 



(156) 



where ip a quark field; m q is the quark mass; tpj^ A B X B ip — A B ipi a i y£gX B j 'ipj0. In Eq. 1156|) we took into account 
that according to Section liV Al 



0, B = l,- 
A*B = 8 



(157) 



Eq. 1|157|) nieans that A B ,B = 1,.-- ,7 are quantized degrees of freedom (disordered phase) but A 8 remains in 
classical state (ordered phase). The field equations for this system are 

"3 



^ 7 2 



2 - 2 ^ Q r + a 2 U m <t> m -(4) 



d v h^ = - A g 2 k\ (rr) b» - {m 2 f v K + g Affix A V, 
[f> (p^ + gAl\ 8 ) ~ m q ] 4, = 0. 

For the simplicity we choose ij)\ a — ip 2a — ^3q- F° r this choice the term rpj^ X 8 ip in Eq. I|160|) vanishes. 
The solution of Eq's (|158fl - 1161|) we search in the form 



(158) 
(159) 

(160) 
(161) 



&„ = ^{/(r), 0,0,0}; ^ = ^(r); 



1 1 



(162) 



\m~ ) = diag{m , 0,0,0} 

here r, 9, ip are spherical coordinates. The corresponding color (but Abelian) radial electric field is 

?8 = _dh 
dr ' 



E 8 



(163) 
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The ansatz for the spinor field is in the standard form 

\ {-l) 1± ^~ L h(r)n j i lm 

here I = j + 1/2, I' = 2j — 1. Thus we have the following equations set 

„ 2 



x 

x" + -x' 

x 



4> [x 2 + Ai 



Mi) ' / 2 



= X [0 2 + A 2 (x 2 - , 
/" + -/ = / 



^1 \ ,2 2 

-) -mo 



1 + x 

a; 

1 - K 



~ (E + m- gf)h = 0, 
h + (E-m-gf)f = 



here we have introduced the dimensionless coordinate x — r<fi(0); redefined v — * v/(f>o, h — ► /i/</>o and 



+ ifj = Z + | 



ifj = / 



Eq's l|165|H167fl nave been solved in Section The numerical calculations of Eq's (|168|l l|l()9|) (with x = 
presented in Fig. ITT1 



(164) 

(165) 
(166) 

(167) 

(168) 
(169) 

(170) 
T) are 



IX. DISCUSSION AND CONCLUSIONS 

In this work we have investigated color defects in a gauge condensate. The gauge condensate is defined as quantized 
non-perturbative degrees of freedom of the SU(3) gauge field (disordered phase). The color defects are residual degrees 
of freedom remaining in a classical state (ordered phase). Using some assumptions and approximations for 2 and 4- 
points Green's functions the SU(3) Lagrangian is reduced to an effective Lagrangian. The 2-point Green's function 
of quantum degrees of freedom is presented as a bilinear combination of scalar fields. The 4-point Green's function is 
presented as a bilinear combination of 2-points Greens functions. The effective Lagrangian describes a gauge theory 
with broken gauge symmetry interacting with scalar fields. The gauge field describe the ordered phase and scalar 
fields - disordered phase. 

In this model we have a few undefined parameters which are connected with a non-linear interaction between initial 
gauge fields. Some of them are defined as eigenvalues at which regular solutions of corresponding field equations 
exist. The obtained solutions describe: infinite flux tube filled with longitudinal color electric and magnetic fields; a 
color electric hedgehog; a defect having two color electric dipoles lying on a line and color magnetic dipoles; two color 
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electric dipoles; two color magnetic dipoles; parton-like defect which can be considered as a bag for the quark held. 
We assume that these defects are quantum excitations in a gauge condensate. 
Finally we would like to list the main results of the paper: 

• An approximate method for non-perturbative quantization of QCD is proposed. 

• The proposed method is based on the decomposition of SU(3) degrees of freedom on ordered and disordered 
phases. 

• The ordered phase is considered as classical gauge helds. This phase forms color defects inside of the disordered 
phase. For example, it can be: a flux tube filled with a longitudinal electric held, a hedgehog filled with a radial 
electric held and so on. 

• The disordered phase is quantum degrees of freedom which squeezes the ordered phase into color defects. 

• Probably the disordered phase is formed by Abelian monopoles. 

• The ordered phase is described by Yang-Mills equations with broken gauge symmetry. 

• The equations for the disordered phase are an analog of Ginzburg - Landau equation for the wave function of 
Cooper pairs. 
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